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A  Predicted  Geomagnetic  Field 
Model  for  Epoch  1990.0 


1.0  Introduction 


2.0  Epoch  1990.0 
Predicted  Field  Model 


The  U.S.  Navy  is  responsible  within  the  Department  of  Defense  for 
production  of  global  geomagnetic  field  models.  These  models  are  used 
for  navigation  and  related  applications,  such  as  Attitude/Heading  Reference 
Systems  (AHRS).  The  Naval  Oceanographic  Office  (NAVOCEANO) 
accumulates  necessary  magnetic  field  measurement  data  and  produces 
such  models  every  5  years.  Because  the  geomagnetic  field  changes  with 
time,  a  model  of  the  time  variation  (secular  variation)  of  the  main 
(core-generated)  field  is  also  produced  for  use  in  updating  the  main  field 
model  in  the  period  between  production  of  field  models. 

The  Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
(NOARL)  is  engaged  in  an  ongoing  basic  research  program  aimed  at 
increasing  understanding  of  time  variations  of  the  geomagnetic  field 
and  determining  the  best  methods  for  modeling  and  predicting  the 
geomagnetic  field.  This  research  has  not  been  completed;  however, 
models  of  secular  variation  have  been  produced  for  each  year  from 
1962.0  through  1983.0  using  vector  magnetic  field  data  at  73  magnetic 
observatories  as  the  data  source.  The  year  is  represented  in  decimal 
form  so  that  1962.0  is  the  start  of  January  1, 1962  in  universal  (Greenwich) 
time. 

An  epoch  1990.0  geomagnetic  field  model  will  be  produced  by 
NAVOCEANO  using  magnetic  field  measurements  from  the  recently 
launched  (April,  1990)  Polar  Orbiting  Geomagnetic  Survey  (POGS) 
satellite  as  part  of  the  input  data.  To  aid  in  the  production  and  evaluation 
of  this  model,  NAVOCEANO  requested  that  NOARL  produce  a  predicted 
field  model  for  1990.0. 

Section  2  of  this  technical  report  describes  production  of  the  predicted 
epoch  1990.0  field  model  from  a  standard  1980.0  field  model  (the 
E)efinitive  Geomagnetic  Reference  Field  model  described  by  Langel,  1988) 
updated  to  1990.0  through  use  of  the  NOARL  secular  variation  modeis. 
In  section  3,  the  NOARL  secular  variation  models  are  described.  An 
error  analysis  for  the  predicted  field  model  is  given  in  section  4. 
Conclusions  given  in  section  5  include  suggestions  on  how  the  predicted 
field  model  might  be  improved. 

A  discussion  of  geomagnetism,  geomagnetic  field  modeling,  and 
International  Geomagnetic  Reference  Fields  is  given  by  Langel  (1987). 


Spherical  harmonic  coefficients  (also  called  Gaussian  coefficients) 
for  the  NOARL  epoch  1990.0  predicted  field  model  are  shown  in  Table  1. 
This  model  is  called  NOARL90.  The  coefficients  are  for  the  Schmidt 
semi-normalized  spherical  harmonic  representation  commonly  used  in 
geomagnetism  by  international  agreement  circa  1939.  The  NOARL90 
model  was  produced  by  updating  the  DGRF  model  for  1980.0  by  use 
of  a  NOARL  secular  variation  model  for  1980.0.  Coefficients  of  the 
DGRF  model  are  shown  in  Table  2.  This  DGRF  model  is  called  DGRF80 
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Table  1.  NOARL90  epoch  1990.0  predicted  field  model. 


N 

G(N,M) 

H(N.M) 

N 

M 

G(N,M) 

H(N,M) 

1 

0 

-29755  9 

0.0 

7 

5 

5.9 

18.0 

1 

1 

-1836.5 

5428.9 

7 

6 

8.1 

-19.9 

2 

0 

-2166  4 

0.0 

7 

7 

-2  5 

-8  3 

2 

1 

3058  1 

-2265  0 

8 

0 

21.0 

0.0 

2 

2 

1720  9 

-434  8 

8 

1 

5.3 

6.4 

3 

0 

1298  8 

0.0 

8 

2 

0  1 

-21.6 

3 

1 

-2240  2 

-298.1 

8 

3 

-11.4 

5.3 

3 

2 

1230.6 

299.7 

8 

4 

-10.4 

-23  8 

3 

3 

825  9 

-326.9 

8 

5 

4.6 

11.2 

4 

0 

931  5 

0.0 

8 

6 

5.8 

11.9 

4 

1 

773  4 

253  5 

8 

7 

4.0 

-16.7 

4 

2 

328  1 

-245  9 

8 

8 

-2.7 

-15.3 

4 

3 

-431.6 

87.5 

9 

0 

5.4 

0.0 

4 

4 

141  9 

-303  2 

9 

1 

8  5 

-20.7 

5 

0 

-214  5 

0  0 

9 

2 

1.0 

16.4 

5 

1 

354  6 

54.7 

9 

3 

-11.6 

10.5 

5 

2 

248  7 

148.3 

9 

4 

8.5 

-6.5 

5 

3 

-110.0 

-160.8 

9 

5 

-4.4 

-7.5 

5 

4 

-168  8 

-73  1 

9 

6 

-12 

6.6 

5 

5 

-45.1 

96.7 

9 

7 

8  2 

8  6 

6 

0 

54  9 

0.0 

9 

8 

3  3 

—6.5 

6 

1 

67.9 

-17  5 

9 

9 

-6  1 

2.0 

6 

2 

61  8 

86.8 

10 

0 

-4  0 

0.0 

6 

3 

-178  7 

66.2 

10 

1 

-4.0 

1.0 

6 

4 

123 

-46.3 

10 

2 

20 

0.0 

6 

5 

15.1 

2.5 

10 

3 

-5.0 

-3.0 

6 

6 

-106  4 

23  5 

10 

4 

-2.0 

6.0 

7 

0 

789 

0.0 

10 

5 

5  0 

-4.0 

7 

1 

-64.0 

-84  4 

10 

6 

3  0 

0.0 

7 

2 

3.7 

-26.4 

10 

7 

1.0 

-10 

7 

3 

24  2 

2  5 

10 

8 

2  0 

4.0 

7 

4 

-1.5 

21.8 

10 

9 

3.0 

0.0 

10 

10 

0.0 

-6  0 

and  is  described  by  Langel  (1988).  The  NOARL  secular  variation  model 
is  called  NOARL80S  and  its  coefficients  are  shown  in  Table  3.  The 
DGRF80  model  was  updated  to  produce  the  predicted  NOARL90  model 
by  adding  coefficients  of  the  NOARL80S  model  multiplied  by  10  years 
to  corresponding  coefficients  of  the  DGRF80  model.  Units  for  the 
NOARL90  and  DGRF80  models  are  nT  while  units  for  the  NOARL80S 
model  are  nT/yr. 

The  NOARL80S  secular  variation  model  was  produced  by  averaging 
NOARL  secular  variation  models  for  years  1977.0  through  1983.0;  that 
is,  corresponding  coefficients  for  these  models  were  averaged.  Since 
the  NOARL  secular  variation  models  for  the  individual  years  1977.0 
through  1983.0  are  in  geomagnetic  coordinates,  it  was  necessary  to 
transform  coefficients  from  the  geomagnetic  coordinate  system  to  the 
standard  geocentric  coordinate  system  aligned  with  Earth’s  rotation  axis. 

Another  predicted  1990.0  field  model  was  produced  from  the 
International  Geomagnetic  Reference  Field  (IGRF)  IGRF85  main  field 
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Table  2.  DGRF80  epoch  1980.0  main  field  model. 


N 

M 

G(N.M) 

H(N.M) 

N 

M 

G(N.M) 

H(N,M) 

1 

0 

-29992 

0 

7 

5 

1 

18 

1 

1 

-1956 

5604 

7 

6 

11 

-23 

2 

0 

-1997 

0 

7 

7 

-2 

-10 

2 

1 

3027 

-2129 

8 

0 

18 

0 

2 

2 

1663 

-200 

8 

1 

6 

7 

3 

0 

1281 

0 

8 

2 

0 

-18 

3 

1 

-2180 

-336 

8 

3 

-11 

4 

3 

2 

1251 

271 

8 

4 

-7 

-22 

3 

3 

833 

-252 

8 

5 

4 

9 

4 

0 

938 

0 

8 

6 

3 

16 

4 

1 

782 

212 

8 

7 

6 

-13 

4 

2 

398 

-257 

8 

8 

-1 

-15 

4 

3 

-419 

53 

9 

0 

5 

0 

4 

4 

199 

-297 

9 

1 

10 

-21 

5 

0 

-218 

0 

9 

2 

1 

16 

5 

1 

357 

46 

9 

3 

-12 

9 

5 

2 

261 

150 

9 

4 

9 

-5 

5 

3 

-74 

-151 

9 

5 

-3 

-6 

5 

4 

-162 

-78 

9 

6 

-1 

9 

5 

c 

-48 

92 

9 

7 

7 

10 

6 

0 

48 

0 

9 

8 

2 

-6 

6 

1 

66 

-15 

9 

9 

-5 

2 

6 

2 

42 

93 

10 

0 

-4 

0 

6 

3 

-192 

71 

10 

1 

-4 

1 

6 

4 

4 

-43 

10 

2 

2 

0 

6 

5 

14 

-2 

10 

3 

-5 

3 

6 

6 

-108 

17 

10 

4 

-2 

6 

H 

0 

72 

0 

10 

5 

5 

-4 

B 

1 

-59 

-82 

10 

6 

3 

0 

B 

2 

2 

-27 

10 

7 

1 

-1 

B 

3 

21 

-5 

10 

8 

2 

4 

B 

4 

-12 

16 

10 

9 

3 

0 

■ 

10 

10 

0 

-6 

model  and  IGRF85S  secular  variation  model  by  linear  extrapolation. 
These  IGRF  models  are  described  in  the  article  by  Langel  (1988).  The 
new  predicted  model  produced  at  NOARL  from  the  IGRF85  and  IGRF85S 
models  was  called  NEW90.  A  difference  model  D90  was  computed 
from  NOARL90  and  NEW90  by  subtracting  corresponding  coefficients 
of  NEW90  from  those  of  NOARL90.  Field  components  for  D90  were 
computed  on  a  grid  10  degrees  by  10  degrees  in  latitude  and  longitude 
and  the  scalar  field  for  D90  was  then  computed  at  the  grid  points  from 
the  square  root  of  the  sum  of  the  squares  of  the  field  components.  The 
scalar  field  for  D90  is  represented  by  the  contour  plot  of  Figure  1  as  a 
function  of  geographical  location;  the  contour  plot  was  derived  from 
the  scalar  field  at  the  grid  points.  Maximum  difference  field  magnitude 
is  about  480  nT  and  occurs  in  the  South  Pacific  at  a  location  far  from 
the  nearest  magnetic  observatory.  Difference  field  D90  reflects  errors 
in  one  or  both  of  the  predicted  field  models.  It  will  be  shown  in  section  4 
that  the  large  difference  field  in  the  South  Pacific  is  probably  due 
mostly  to  errors  in  the  NEW90  predicted  model. 
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Table  3.  NOARL80S  epoch  1980.0  secular  variation  model. 


3.0  NOARL  Secular 
Variation  Model 


N 

M 

G(N.M) 

H(N,M) 

N 

M 

G(N,M) 

1 

0 

23  608 

0.000 

7 

0 

0  690 

0  000 

1 

1 

1 1  954 

-17.508 

7 

1 

-0.500 

-0.237 

2 

0 

-16  939 

0.000 

7 

2 

0.173 

0  057 

2 

1 

3  113 

-13  599 

7 

3 

0.320 

0.749 

2 

2 

5  787 

-23  475 

7 

4 

1.048 

0.578 

3 

0 

1  776 

0.000 

7 

5 

0.492 

0  004 

3 

1 

-6  019 

3  791 

7 

6 

-0.286 

0.312 

3 

2 

-2.037 

2  874 

7 

7 

-0.055 

0  172 

3 

3 

-0  709 

-7  492 

8 

0 

0  298 

0.000 

4 

0 

-0.653 

0.000 

8 

1 

-0.067 

-0.059 

4 

1 

-0.860 

4.153 

8 

2 

0.005 

-0.356 

4 

2 

-6  992 

1.111 

8 

3 

-0.043 

0.133 

4 

3 

-1  256 

3  452 

8 

4 

-0.339 

-0.179 

4 

4 

-5.709 

-0.616 

8 

5 

0  056 

0.221 

5 

0 

0  350 

0  000 

8 

6 

0  282 

-0.415 

5 

-0.241 

0.868 

8 

7 

-0  201 

-0.371 

-1  231 

-0.168 

8 

8 

-0  174 

-0.033 

3 

-3  595 

-0.977 

9 

0 

0.045 

0.000 

5 

4 

-0.681 

0.490 

9 

1 

-0.148 

0  032 

5 

5 

0  292 

0.474 

9 

2 

0.002 

0.036 

6 

0 

0.688 

0.000 

9 

3 

0  045 

0.147 

6 

1 

0.189 

-0.245 

9 

4 

-0.049 

-0.152 

6 

2 

1.978 

-0.624 

9 

5 

-0.138 

-0.148 

6 

3 

1.330 

-0  483 

9 

6 

-0.020 

-0  245 

6 

4 

0.825 

-0.325 

9 

7 

0.120 

-0.141 

6 

5 

0.111 

0.454 

9 

8 

0  128 

-0.053 

6 

6 

0.157 

0.651 

9 

9 

-0.107 

-0.002 

A  set  of  secular  variation  models  has  been  produced  by  NOARL 
from  annual  means  of  magnetic  observatory  vector  field  measurements. 
The  models  were  produced  from  first  differences  of  annual  means  of 
field  components  at  73  observatories  which  had  vector  measurements 
for  each  of  the  years  1961.5  through  1983.5.  A  separate  secular  variation 
model  was  produced  for  each  year  1962.0  through  1983.0.  The  models 
were  of  spherical  harmonic  degree  and  order  9  for  internal  source  terms 
and  the  Gaussian  coefficients  were  obtained  for  a  geomagnetic  coordinate 
system.  External  source  terms  were  included  in  the  models  for  zonal 
spherical  harmonics  of  degrees  1.  3,  5,  7. 

The  method  of  stochastic  inversion  was  used  to  produce  the  models; 
this  method  is  discussed  by  McLeod  (1986).  Use  of  this  method  requires 
an  error  covariance  matrix  for  the  first  differences  of  the  observatory 
field  components  and  an  a  priori  parameter  covariance  matrix.  Both  of 
these  covariance  matrices  were  taken  to  be  diagonal,  corresponding  to 
an  assumed  statistical  independence  for  the  different  measurement  errors 
and  assumed  isotropy  and  spatial  stationarity  for  the  random  process 
rcprc.seniing  the  internal  source  geomagnetic  field.  Parameter  variances 
for  the  internal  source  field  were  based  on  a  function  that  was  fitted  by 
McLeod  (1985)  to  the  1975  secular  variation  spatial  power  spectrum 
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Figure  1.  Magnitude  of  difference  field  at  1990.0  for  predicted  field  models. 
NOARL90  -  NEW90.  Units  are  nT. 

derived  from  the  GSFC(9/80)  model  described  by  Langel  et  al.  (1982). 
The  function  fit  to  the  1975  secular  variation  spatial  power  spectrum  was: 

Fin)  =  1 00  (2n  +  1 9  ida)^  (nT/yr)‘  ( 1 ) 

where  n  is  degree  of  the  spherical  harmonic,  c  is  core  radius  (3485  km), 
a  is  mean  Earth  radius  (6371.2  km).  This  function  Fin)  represents  the 
mean  squared  value  of  that  portion  of  the  field  derivative  represented 
by  spherical  harmonics  of  degree  n,  and  was  given  by  McLeod  (1985). 
Since  there  are  i2n  +  1)  spherical  harmonic  coefficients  of  degree  n  and 
the  mean  squared  field  associated  with  each  coefficient  is  (n+  1) 
multiplied  by  the  square  of  the  coefficient,  the  a  priori  parameter 
variances  used  in  the  stochastic  inversion  were: 

Cin)  =  400  in  +  l/2)\n  +  l)-'(c/a)2"  (nT/yr)^  (2) 

where  Cin)  is  the  a  priori  variance  for  any  of  the  coefficients  of  degree  n. 

Locations  of  the  73  magnetic  observatories  used  for  the  secular  variation 
models  are  shown  in  Figure  2a.  Table  4  lists  estimated  rms  errors  for 
each  field  component  at  each  observatory.  These  error  estimates  are 
based  on  rms  values  of  third  differences  of  individual  field  components 
at  individual  observatories  for  the  interval  1961.5-1983.5.  The  field 
components  arc  in  geomagnetic  coordinates  and  a  model  external  field  was 
subtracted  from  observatory  measurements  before  computing  third 
differences.  A  more  detailed  discussion  of  the  production  of  the  NOARL 
secular  variation  models  is  given  in  a  paper  by  McLeod  (in  preparation). 
The  a  priori  variances  used  in  the  inversion  for  the  external  zonal 
coefficients  of  degrees  1,  3,  5,  7  were  178,  0.44,  0.1 1, 0.1 1,  respectively 
in  units  of  (nT/yr)^. 

Table  4  shows  estimated  rms  errors  in  the  model  at  various  observatories 
in  units  of  nT/yr  and  also  expressed  as  a  percentage  of  measurement 
error  at  the  various  observatories.  Estimated  model  error  was  computed 
from  the  error  covariance  matrix  and  a  priori  parameter  covariance 


A  Predicted  Geomagnetic  Field  Model  for  Epoch  1990.0 


5 


Figure  2.  (a)  Locations  of  magnetic 
observatories  used  for  NOARL80S 
secular  variation  model  and 
(b)  estimated  magnitude  of  vector  field 
error  for  NOARL  secular  variation 
models.  Seven  of  these  secular 
variation  models  were  averaged 
to  produce  NOARL80S.  Units  are 
nTlyr. 


4.0  Error  Analysis 


matrix  by  the  method  discussed  by  McLeod  (1986).  Estimated  model 
error  was  also  computed  on  a  grid  of  latitude  and  longitude.  Estimated 
rms  error  of  the  model  vector  field  is  shown  as  a  function  of  geographic 
location  in  the  contour  plot  of  Figure  2b. 

It  can  be  seen  from  the  contour  plot  of  Figure  2b  and  the  estimated 
errors  of  Table  4  that  the  model  is  most  accurate  in  regions  where  there 
is  a  high  density  of  observatories,  such  as  Europe.  The  model  is  most 
inaccurate  in  regions  far  from  a  magnetic  observatory,  such  as  the  South 
Pacific.  At  isolated  observatories,  such  as  Honolulu  or  Hermanus,  the 
model  closely  fits  measurement  data,  so  the  model  error  is  nearly  the 
same  as  measurement  error;  but  in  regions  where  there  are  many 
observatories,  model  error  is  less  than  observatory  measurement  error 
because  the  model  effectively  averages  measurements  from  several  nearby 
ob.se  rvaiories. 


Errors  in  the  NOARL90  predicted  field  model  are  due  to  three  causes: 
(a)  errors  in  the  DGRF80  model;  (b)  errors  in  the  NOARL80S  secular 
variation  model;  and  (c)  errors  due  to  the  fact  that  secular  variation  is 
not  constant  from  1980  to  1990  but  is  time  varying.  Errors  due  to  errors 
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Table  4.  Estimated  rms  measurement  errors  and  model  errors  at  73  magnetic  observatories 
Units  are  nT/vr.  Letters  X.  Y,  and  Z  refer  to  north,  east,  and  vertical  field  components,  respectively. 
Field  components  are  in  geomagnetic  coordinates. 


Observatory 

Geographic 

Lat.  Long. 

Meas.  Error 

X  Y  Z 

Model  Error 

X  Y  Z 

Percent 

X 

Model 

Y 

Error 

Z 

Alibag 

16.638 

72.872 

15.8 

5  2 

19.7 

5  0 

4.1 

7  3 

31.6 

78.2 

39  6 

Almeria 

36.853 

-2.460 

4.1 

2  2 

5  4 

1.8 

1.4 

2  3 

44.6 

63.5 

42.1 

Apia 

-13.807 

-171.775 

2.7 

2  9 

4.5 

2  6 

2.9 

4  4 

96.7 

97.5 

98,0 

Argentine  Island 

-65.245 

-64  258 

3.5 

2.8 

2  9 

3.3 

2  8 

2  8 

94  8 

98.1 

39  3 

Baker  Lake 

64.333 

-96.033 

126 

29  7 

14  6 

2.2 

2.1 

3  3 

17.3 

6.9 

22  9 

Bangui 

4.437 

18.565 

4.6 

9.4 

7  4 

3  8 

6.1 

6.4 

82.0 

65.6 

86.2 

Beijing 

40,040 

116  175 

4.1 

1.8 

6  0 

1.6 

1.4 

2,4 

40  0 

77.4 

40.3 

Cape  Wellen 

66.163 

-169.835 

8.8 

13  4 

8  9 

2  9 

2.4 

2.5 

33.5 

17.8 

28  " 

Cha  Pa 

22.350 

103  833 

15.0 

9  7 

24.8 

4.7 

3  8 

5.9 

31  4 

39.4 

23.6 

Chambon-La-Foret 

48.023 

2  260 

2  3 

2  4 

3.2 

0  6 

0.6 

0  7 

24.8 

24.1 

21.8 

College 

64  860 

-147.837 

3.2 

1.2 

2  6 

2  0 

11 

2.0 

63.9 

91.0 

75.3 

Dixon  Island 

73.543 

80.562 

15.6 

5  4 

16.7 

3.0 

2  9 

4  1 

19.3 

53.7 

24  3 

Dombas 

62.073 

9  117 

2  1 

2  1 

1.8 

0.7 

0.5 

0  8 

32.9 

23  8 

41.9 

Dourbes 

50.097 

4  595 

11 

1.7 

2.3 

0  4 

0  5 

0.6 

36  9 

27  1 

24  5 

Dusheti 

42  092 

44  705 

15 

1.5 

2  2 

1.3 

1.4 

18 

91  0 

90.7 

83  2 

Eskdalemuir 

55.317 

-3.200 

1.6 

12 

6  3 

0  5 

0  5 

0  7 

33  9 

41.7 

11.5 

Fredericksburg 

38.205 

-77.373 

2  2 

12 

2.6 

1,9 

1.2 

2  5 

86  7 

96  2 

95  8 

Furstenfeldbruck 

48  165 

11.277 

2  2 

1.1 

1.0 

0.4 

0.4 

0  5 

18  1 

38.4 

54  3 

Gnangara 

-31.783 

1 1 5 .950 

2  0 

15 

2  2 

2  0 

1.5 

2  2 

97.9 

98  8 

99  4 

Godhavn 

69  252 

-53.533 

2.6 

1.6 

6  7 

1.9 

1.3 

3  7 

73  1 

83  7 

55  2 

Gonotayezhnaya 

43.683 

132.167 

2.5 

29 

11.9 

1.1 

1.1 

1.8 

46.4 

38  1 

14.7 

Grocka 

44.633 

20.767 

1.3 

1.0 

2.4 

0  6 

0.5 

0.9 

44  0 

51.2 

35  9 

Guam 

13.583 

144.870 

5.7 

1.8 

3.3 

3.4 

1.7 

2.9 

60  1 

93.3 

88  5 

Hartland 

50  995 

-4.483 

1.7 

2  1 

1.5 

0  6 

0  6 

0  8 

35.1 

28.3 

53.3 

Hel 

54  608 

18.815 

2.5 

2.7 

2  6 

0.4 

0.4 

0  6 

16  3 

14.9 

24  1 

Hermanus 

-34.425 

19.225 

12 

1.0 

4.0 

12 

1.0 

3.7 

98  8 

99  3 

91  9 

Honolulu 

21.320 

-158.002 

2.7 

14 

2  4 

2.6 

14 

2.3 

94  3 

99  1 

99  2 

Huancayo 

-12  045 

-75.340 

2  6 

25 

3.1 

24 

2  4 

3,0 

92.1 

969 

97.3 

Hurbanovo 

47.873 

18.190 

2  6 

1.8 

3.0 

0.4 

0.4 

0.6 

16.6 

22  9 

20  3 

Kakioka 

36.230 

140  190 

0.4 

0  6 

0.8 

0  4 

0.5 

0  6 

92  7 

86.5 

70.0 

Kanoya 

31  420 

130  882 

1.9 

13 

0  9 

1.0 

0.9 

0.3 

52  2 

69.6 

87  0 

Kanozan 

35,253 

139.960 

1.4 

1.6 

0  9 

0.4 

0,5 

0.6 

28  5 

320 

64  4 

L'  Aquila 

42.383 

13.317 

1.6 

1.8 

1.5 

0  8 

0.7 

0.9 

49.2 

40.7 

57.1 

Leirvogur 

64.183 

-21  .700 

4.1 

1.7 

2.7 

1.5 

1.1 

1.7 

36.2 

66. 1 

62.5 

Lerwick 

60  133 

-1 .183 

1.2 

18 

1.5 

0.7 

0.6 

0  8 

61.1 

30.7 

55.1 

Lovo 

59  345 

17.827 

1.2 

0.6 

1.0 

0  5 

0  4 

0.6 

46  8 

74.3 

64.5 

Lvov 

49.900 

23.750 

6  0 

2  1 

7  2 

0.5 

0.5 

0  8 

8.2 

24.0 

10.6 

M'  Bour 

14.392 

-16.958 

3  5 

3  4 

3  3 

3  2 

3,1 

3.1 

91.3 

91.5 

95  5 

Macquarie  Island 

-54.500 

158  950 

3  3 

8  0 

7.4 

3.1 

6.8 

7.0 

95  4 

84  5 

95.0 

Memambetsu 

43  907 

144  193 

2.0 

1.1 

2.9 

1.3 

0.9 

14 

62,3 

81.1 

48.6 

Mirny 

-66  550 

93.017 

8  4 

11.3 

19.3 

4.6 

4.7 

8  3 

54  7 

42.1 

42  8 

Muntinlupa 

14.375 

121.015 

6  8 

7  7 

12  3 

4  1 

3,7 

5.1 

59.6 

47.5 

41.3 

Nagycenk 

47.633 

16  717 

8  7 

5  5 

7.0 

0  4 

0  4 

0.6 

4,9 

7  5 

8  5 

Niemegk 

52  072 

12.675 

0  4 

1.0 

3.1 

0  3 

0.4 

0  5 

76.2 

34  3 

170 

Nurmijarvi 

60  508 

24.655 

9.7 

36 

6.9 

0  7 

0.6 

0.8 

7,1 

17.2 

11,3 

Panagyurishte 

42  515 

24  177 

1.2 

14 

3.3 

0  8 

0.6 

1.1 

69.3 

45  3 

34.2 

Patrony 

52.167 

104  450 

2  0 

3  6 

3  9 

18 

2  6 

3.3 

86  3 

70  6 

83  5 

Pilar 

-31.667 

-63.883 

10.9 

3  1 

4.3 

5  0 

2  8 

4.0 

45.6 

88.9 

93.0 

Pleshenitzi 

54.500 

27.883 

3  9 

7  4 

4  0 

0.6 

0  7 

0  9 

14.2 

9  5 

22  8 
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Table  4  (cent.) 


Observatory 

Geographic 

Lat.  Long 

Meas.  Error 

X  Y  2 

Model  Error 

X  Y  2 

Percent 

X 

Model 

Y 

Error 

2 

Port  Moresby 

-9.408 

147.152 

6.1 

1.5 

8.3 

4.5 

1.4 

7.0 

73  3 

98.0 

84.8 

Pcrt-Aux-Francais 

-49.350 

70.200 

5.5 

4  6 

1.7 

4  6 

3.8 

1.6 

83.6 

83.1 

99  0 

Resolute  Bay 

74.700 

-94.900 

2.5 

4  7 

7.5 

1.7 

2.4 

4  4 

70  0 

50.5 

59  3 

San  Juan 

18.113 

-66.150 

2.4 

1.8 

3  0 

2  2 

1.7 

2.8 

91  1 

97.0 

94  1 

Sheshan 

3  (.097 

121.187 

2  3 

2  6 

4.7 

1.6 

1.4 

18 

66  8 

55.9 

37  8 

Sitka 

57.058 

-135.325 

2.7 

0  9 

3.0 

1.7 

0.8 

1.7 

62.9 

89.1 

56.7 

Sodankyla 

67.368 

26.630 

2  4 

12 

2.0 

1.3 

0.9 

1.2 

53.1 

80  3 

60  7 

Stepanovka 

46.783 

30.883 

2.4 

1.8 

3  4 

0.7 

0  8 

1.1 

28  4 

41.8 

31.1 

Surlari 

44  680 

26  253 

12  3 

3.8 

20  0 

0.7 

0  6 

1.0 

5  5 

16.4 

5.1 

Tamanrasset 

22.792 

5.527 

7.5 

20.5 

16.5 

3.3 

2.7 

4  8 

44.3 

13.1 

28  9 

Tihany 

46.900 

17.893 

6  1 

4.3 

20.0 

0.5 

0  4 

0.6 

7.6 

10.0 

3.2 

Trivandrum 

8.483 

76  950 

6.7 

10.7 

8.8 

5.6 

5  8 

7.5 

83  1 

53.8 

85.3 

Tucson 

32.247 

-110  833 

4  1 

3  1 

4.7 

3.6 

2.6 

3.8 

86  9 

85  1 

81.7 

Valentia 

51  933 

-10  250 

1.4 

0  9 

19 

0  8 

0.8 

1.0 

55.0 

81.9 

52  3 

Vassouras 

-22.400 

-43  650 

124 

7.0 

5  0 

6  8 

5.2 

4.6 

55  1 

74  9 

91.2 

Victoria 

48.517 

-123  417 

2.6 

16 

2  6 

18 

13 

19 

70.5 

81.9 

72.9 

Voyeykovo 

59.950 

30  705 

3  9 

3  0 

3  9 

0  8 

0.9 

10 

20  1 

29  5 

25  4 

Wien-KobenzI 

48  265 

16.318 

1.0 

1.4 

14 

0.4 

0  4 

0  6 

42  4 

29  3 

41  .5 

Wingst 

53.743 

9.073 

1.5 

0  9 

1.8 

0  3 

0  4 

0.5 

23  2 

38.5 

29  3 

Witteveen 

52.813 

6.668 

6.0 

17 

1.3 

0.4 

0.4 

0.5 

5.9 

22  7 

42  5 

Yuzhno  Sakhalin 

46  950 

142  717 

6.0 

9.6 

3  4 

1.6 

1.2 

1.8 

27  2 

12.3 

52  1 

in  NOARL80S  can  be  estimated  conservatively  from  the  contour  plot 
shown  in  Figure  2b  by  multiplying  errors  in  nT/yr  shown  in  that  figure 
by  the  10-year  interval  from  1980  to  1990.  This  estimate  is  conservative 
because  errors  shown  in  Figure  2b  are  for  a  model  based  on  first  differences 
of  annual  means  for  a  single  pair  of  years  while  NOARL80S  is  the 
average  of  seven  such  models. 

In  order  to  estimate  accuracy  of  the  NOARL90  model,  a  similar 
model  was  computed  for  epoch  1980.0  and  named  NOARL80.  This  model 
is  an  update  of  the  DGRF70  field  model  computed  with  the  aid  of  a 
secular  variation  model  NOARL70S.  This  secular  variation  model  was 
computed  by  averaging  NOARL  secular  variation  models  for  years  1967.0 
through  1973.0.  Thus,  expected  rms  errors  for  the  NOARL80  model  are 
the  same  as  those  for  the  NOARL90  model,  provided  the  underlying 
random  process  is  temporally  stationary.  (The  term  “expected  rms  errors” 
is  used  in  the  statistical  sense.)  Coefficients  of  the  DGRF80  model  were 
subtracted  from  corresponding  coefficients  of  the  NOARL80  model  to 
form  a  difference  model  D80.  If  errors  in  the  DGRF80  model  are  negligible 
in  comparison  with  errors  in  the  NOARL80  predicted  model,  then  D80 
represents  the  error  for  the  NOARL80  model.  Contour  plots  of  the  field 
derived  from  the  difference  model  D80  are  shown  in  Figures  3a-d. 
These  figures  show  the  scalar  field,  north,  east,  and  vertical  components, 
respectively. 

The  contour  plot  of  Figure  3a  shows  no  similarity  to  the  contour  plot 
of  Figure  2b;  thus,  errors  in  the  NOARL70S  secular  variation  model  are 
not  the  principal  source  of  errors  in  the  NOARL80  predicted  field  model. 
Since  it  is  improbable  that  errors  in  the  DGRF7()  or  DGRF80  models  are 
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Figure  3.  (a)  magnitude,  (h)  north 
component,  (c)  east  component,  and 
id)  vertical  component  of  difference 
field  at  1980.0  for  predicted  field 
model  NOARL80  minus  standard 
model  DGRF80.  Units  are  nT. 
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large  enough  lo  produce  the  difference  fields  shown  in  Figure  3a,  the 
difference  field  D80  must  be  due  to  changes  in  secular  variation  between 
1970  and  1980.  To  verify  that  this  is  the  case,  coefficients  of  secular 
variation  model  NOARL70S  were  subtracted  from  corresponding 
coefficients  of  secular  variation  model  NOARL80S  to  obtain  change  in 
secular  variation  for  this  lO-year  time  period.  Field  components  of  the 
difference  model  are  shown  in  contour  plots  of  Figures  4a-d.  Contour 
plots  of  Figure  4  arc  very  similar  to  contour  plots  of  Figure  3.  Numerical 
values  for  Figure  3  arc  about  five  times  those  for  Figure  4,  as  would  be 
expected  if  errors  shown  in  Figure  3  arc  due  to  the  change  of  secular 
variation  shown  in  Figure  4, 

Because  NOARL90  and  NOARL80  were  produced  by  the  same 
methods,  most  of  the  error  in  the  NOARL90  predicted  model  is  probably 
due  to  changes  in  secular  variation  during  the  time  period  1980  through 
1990,  and  estimated  rms  error  in  the  field  corresponding  to  the  NOARL90 
model  is  about  the  same  as  rms  errors  seen  in  Figure  3.  Thus,  estimated 
rms  error  for  the  predicted  model  NOARL90  is  200  nT  and  peak  errors 
for  the  field  components  arc  estimated  to  be  +350,  +250,  +500  nT  for 
north,  cast,  and  vertical  components,  respectively.  Estimated  peak  value 
of  the  vector  error  magnitude  is  500  nT. 

Because  the  error  for  the  NOARL80  predicted  model  does  not  show 
a  large  peak  in  the  South  Pacific  in  the  contour  plot  of  Figure  3a,  the 
peak  in  the  South  Pacific  region  seen  in  the  contour  plot  of  Figure  1  is 
probably  not  due  to  errors  in  the  NOARL90  model,  but  instead  the  peak 
is  probably  due  to  errors  in  the  NEW90  model.  Such  errors  would  be 
expected  if  the  spherical  harmonic  degree  and  order  for  the  IGRF85S 
model,  used  to  update  the  IGRF85  model  to  1990.0,  were  too  high  for 
the  density  of  observatories  in  the  South  Pacific  region.  Spherical  harmonic 
models  produced  by  the  usual  least  squares  fit  are  known  to  attain  large 
error  magnitudes  in  regions  where  the  density  of  observatories  is  low  if 
the  degree  and  order  of  the  model  is  too  great.  The  method  of  stochastic 
inversion  u.scd  to  produce  the  NOARL  secular  variation  models  minimizes 
this  type  of  error. 

In  order  to  see  what  rms  errors  would  result  if  the  DGRF80  model 
were  used  to  represent  the  epoch  1990.0  field  without  using  a  secular 
variation  model  for  updating  the  DGRF80  model,  the  coefficients  of  the 
DGRF70  model  were  subtracted  from  conesponding  coefficients  of 
the  DGRF80  model.  The  scalar  field  of  the  difference  model  is  represented 
in  Figure  5  by  a  contour  plot  as  a  function  of  geographical  location.  It 
can  be  seen  by  comparing  Figure  5  with  Figure  3a  that  errors  are  about 
four  times  as  large  as  they  are  when  a  secular  variation  model  is  used 
for  updating  an  earlier  field  model. 

As  previously  noted,  errors  in  the  predicted  field  models  NOARL90 
and  NOARL8()  arc  mainly  due  to  changes  in  secular  variation  in  the 
10-ycar  period  between  the  epoch  of  the  predicted  model  and  the  epoch 
of  the  model  being  updated.  Thus,  it  was  thought  that  these  errors  might 
be  reduced  if  the  secular  variation  model  were  first  updated  to  the 
midpoint  of  the  lO-ycar  interval.  A  secular  variation  model  NOARL75S 
was  produced  from  NOARL70S  by  averaging  the  NOARL  secular  variation 
models  for  1967.0  through  1969.0  and  1971.0  through  1973.0,  subtracting 
the  first  average  from  the  second  and  dividing  by  4  years  to  get  the 
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Figure  4.  (a)  magnitude,  (b)  north 
component,  (c)  east  component,  and 
(d)  vertical  component  of  lO-year 
accumulated  difference  field  for 
secular  variation  model  NOARL80 
minus  secular  variation  model 
NOARL70S.  Units  are  nT/yr. 
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Figure  5.  Change  of  geomagnetic 
field  1970-1980.  Magnitude  of 
difference  field  for  model  DGRF80 
minus  model  DGRF70.  Units  are  nT. 


5.0  Summary  and 
Conclusions 


secular  acceleration  for  1970.0,  multiplying  the  secular  acceleration  by 
5.0  to  get  the  change  of  secular  variation  from  1970.0  to  1975.0,  and 
finally  adding  coefficients  for  the  change  of  secular  variation  to 
corresponding  coefficients  of  NOARL70S  to  get  NOARL75S.  The  new 
secular  variation  model,  NOARL75S,  was  u.scd  to  update  the  DGRF70  ^ 

model  to  1980.0  in  the  same  way  that  NOARL70S  was  previously  used 
to  obtain  the  predicted  model  NOARL80.  The  new  predicted  model  was 
called  NOARL80-2.  The  difference  D80-2  between  the  predicted  model 
NOARL80-2  and  the  DGRF80  model  was  found.  The  rms  field 
corresponding  to  D80-2  was  larger  than  the  rms  field  corresponding  to 
D80;  therefore,  this  new  method  for  constructing  a  predicted  field  model  # 

is  apparently  not  as  good  as  the  method  used  to  produce  NOARL80  and 
NOARL90. 


A  predicted  geomagnetic  field  model  for  epoch  1990.0  has  been 
developed.  This  model  is  called  NOARL90  and  is  of  spherical  harmonic  ^ 

degree  and  order  10.  Model  NOARL90  is  based  on  the  DGRF80  model 
updated  to  1990.0  through  use  of  a  secular  variation  model  NOARL80S 
for  epoch  1980.0.  Estimated  rms  error  for  NOARL90  is  200  nT  while 
peak  errors  are  estimated  to  be  approximately  500  nT  in  magnitude  for 
the  vector  field  and  +350,  +250,  +500  nT  for  the  north,  east,  and  vertical  • 

field  components,  respectively.  Errors  are  primarily  due  to  changes  in 
secular  variation  during  the  10  years  between  the  epoch  of  secular  variation 
model  NOARL80S  and  the  epoch  of  predicted  model  NOARL90. 

Estimated  errors  would  be  about  four  times  larger  if  DGRF80  were 
u.sed  without  updating  to  represent  the  epoch  1990.0  field. 

A  different  updating  method  that  attempted  to  take  secular  acceleration  * 

into  account  was  investigated;  however,  the  particular  method  tried  resulted 
in  larger  estimated  errors.  Research  to  determine  the  best  method  for 
predicting  the  geomagnetic  field  is  continuing. 

The  secular  variation  model  NOARL80S  was  based  on  secular  variation 
models  produced  at  NOARL  based  on  annual  mean  data  from  73  magnetic  • 

observatories  that  had  continuous  vector  data  from  1961.5  through  1983.5. 

No  observatory  data  later  than  1983.5  were  used  for  any  of  the  models. 
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6.0  References 


There  are  some  large  regions  of  the  Earth  where  there  were  no 
observatories  that  had  continuous  vector  data  for  the  chosen  lime  period. 

More  recent  observ-atory  measurements  are  now  available  and  have 
been  supplied  to  NOARL  by  NAVOCEANO  on  magnetic  tape.  An 
improved  geomagnetic  field  model  could  be  produced  for  epoch  1990.0 
and  model  accuracy  estimated  by  using  much  the  same  approach  that 
was  used  for  the  production  of  NOARL90,  but  using  more  recent 
observatory  data  and  a  larger  number  of  observatories  belter  distributed 
about  the  Earth.  It  is  possible  to  select  a  larger  number  of  observatories  by 
removing  the  requirement  that  the  observatories  have  continuous  vector 
data  back  to  1961.5.  The  requirement  that  the  data  be  continuous  back 
to  1961.5  relates  to  the  NOARL  research  program;  however,  this 
requirement  is  not  necessary  or  desirable  for  the  specific  goal  of  producing 
an  epoch  1990.0  magnetic  field  model. 
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